Universal decay rule for reduced widths 
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Emission processes including a-decay, heavy cluster decays, proton and di-proton emission are 
analyzed in terms of the well known factorisation between the penetrability and reduced width. By 
using a shifted harmonic oscilator plus Coulomb cluster-daughter interaction it is possible to derive a 
linear relation between the logarithm of the reduced width squared and the fragmentation potential, 
denned as the difference between the Coulomb barrier and Q- value. This relation is fulfilled with 
a good accuracy for transitions between ground states, as well as for most a-decays to low lying 
2 + excited states. The well known Viola-Seaborg rule, connecting half lives with the Coulomb 
parameter and the product between fragment charge numbers, as well as the Blendowke scalling 
rule connecting the spectroscopic factor with the mass number of the emitted cluster, can be easily 
understood in terms of the fragmentation potential. It is shown that the recently evidenced two 
regions in the dependence of reduced proton half-lives versus the Coulomb parameter are directly 
connected with the corresponding regions of the fragmentation potential. 

PACS numbers: 21.10.Tg, 23.50.+Z, 23.60.+e, 23.70.+j, 25.70.Ef 
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The family of emission processes triggered by the 
strong interaction contains various decays, namely par- 
ticle (proton or neutron) emission, two-proton emission, 
a-decay, heavy cluster emission and binary or ternary 
fission. There are also other nuclear decay processes in- 
duced by electromagnetic (7-decay) or weak forces 
decays). The purpose of this work is to investigate 
only the first type of fragmentation, where the emitted 
fragments are left in ground or low-lying excited states. 
They are called cold emission processes and are presently 
among important tools to study nuclei far from the sta- 
bility line. Nuclei close to the proton drip line are inves- 
tigated through proton emission, while the neutron drip 
line region is probed by cold fission processes. On the 
other hand superheavy nuclei are exclusively detected 
by a-decay chains Actually the first paper in theo- 
retical nuclear physics applying quantum mechanics Q 
was devoted to the description of the a-decay in terms 
of the penetration of a preformed particle through the 
Coulomb barrier. 

There are two goals of this paper. The first one is to 
explain the well known Viola-Seaborg rule [H, valid for 
all kinds of cold emission processes. It turns out that it 
is possible to give a simple interpretation of this rule in 
terms of two physical quantities, namely the Coulomb pa- 
rameter, connected with the penetrability, and the frag- 
mentation potential, connected with the reduced width. 
An universal linear dependence between the logarithm 
of the reduced width and fragmentation potential will 
be derived. It will be shown that this interpretation is 
valid not only for transitions between ground states, but 
also for transitions to excited states. On the other hand, 
the scalling dependence of spectroscopic factors in heavy 
cluster decays versus the mass numbers of the emitted 
cluster can be also understood in terms of the fragmen- 
tation potential. 



Let us consider a binary emission process P — > D + 
C from a parent (P) to the daughter nucleus (D) and 
the lighter cluster (C), which can be in particular an a- 
particle or a proton. The total decay width is the sum of 
partial decay widths corresponding to different angular 
momenta, given by [J] 



r, = 2fl(£»,r)7?C9,r) , 



(2.1) 



where it was introduced the standard penetrability and 
reduced width squared Q 



Pi{E h r) = 
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The outgoing spherical Coulomb-Hankel function 
Hi (xi-> K l r ) depends upon two variables, namely the 
Coulomb (or twice the Sommerfeld) parameter 



XI = 



Z D Z c e 2 



hvi 



(2.3) 



and the reduced channel radius pi — Kir. Hereu; = hKi/fi 
and Ki = y/2fiEi/h are the asymptotic relative velocity 
and momentum between the emitted fragments, respec- 
tivelly, in terms of the reduced mass of the daughter- 
cluster system p. It is also defined the center of mass 
(cm) channel energy E\ = Q — E\ ex ' of emitted frag- 
ments, in terms of the difference between the total energy 
(Q-value) and the excitation energy of the daughter nu- 
cleus E\ ex \ The internal component s\ mt \/3, r) at a cer- 
tain radius r inside the Coulomb barrier, is for deformed 
emitters a superposition of different Nilsson components 



2 



multiplied by the propagator matrix [||, depending on 
deformation parameters (3, i.e. 



(2.4) 



For spherical emitters with Kiy = Su> it coincides with 
the wave function component fj- int ' (r) . 

The half life is defined by the inverse of the total decay 
width, i.e. 



T = 



Mn 2 



(2.5) 




FIG. 1: Logarithm of half lives for a-decays from even-even 
nuclei versus Coulomb parameter 12. Sj) . Differents lines con- 
nect decays from nuclei with the same charge number. 



Inside the Coulomb barrier the complex Coulomb- 
Hankel function practically coincides with the real irregu- 
lar Coulomb function and has a very simple WKB ansatz 



exp 



X (arccos^fx — \J x(l — x) 



,V4 



H ( Q + \x,p)Ci 



(2.6) 



where, with the external turning point r e = Zi^e 2 /E 
and barrier energy V a = Z 1 Z 2 e 2 /r, there were introduced 
the following notations 



p r 
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(2.7) 



There were considered here for simplicity transitions be- 
tween states with the same angular momentum I, as for 
instance a-decays or proton emission processes between 
ground states. Thus, the logarithm of the so-called re- 
duced half-life, corrected by the exponential centrifugal 
factor squared Cf, defined by the second line of this re- 
lation, i.e. 

T 

logioT re d = logio-^2 
°; 

Mn 2 o 
= lo gio —gr - lo 9w 2p i - l °9wli ,(2.8) 

should be proportional with the Coulomb parameter, i.e. 



logioT re d = a \ + &o 



(2.9) 



Notice that in most of decay processes between ground 
states one has boson fragments, with zero angular mo- 
mentum, i.e. Ci = 1. The case with I ^ is connected 
with fermions, i.e. proton emission, where C; ^ 1. 
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FIG. 2: Logarithm of half lives for heavy cluster decays and the 
corresponding a-decays from the same mother nuclei versus 
Coulomb parameter S2.3}) . Differents symbols denote charge 
number of the emitted cluster. 



The above relation is also called Geiger-Nuttall law, 
discovered in 1911 for a-decays between ground states 
(where the angular momentum carried by the a-particle 
is I = 0). The explanation of this law was given by G. 
Gamow in 1928 [2], in terms of the quantum-mechanical 
penetration of the Coulomb barrier, i.e. the first line of 
Eq. (|2.2[) . It is characterized by the Coulomb parameter, 
which is proportional with the ratio Zjj/ \JQ a . 

The a-decays between ground states are characterized 
by a remarkable regularity, especially for transitions be- 
tween ground states of even-even nuclei. The fact that 
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a-transitions along various isotopic chains lie on separate 
lines, as stated by the Viola-Seaborg rule i.e. 



log w T = 



d + 0.2 



+ hZ L 



(2.10) 



is connected with different a-particle reduced widths, 
multiplying the penetrability in (|2.1|) . From Eq. (|2.8|) 
it becomes clear that the reduced width should depend 
upon the charge of the daughter nucleus. This feature is 
shown in figure [T] Still in doing systematics along neu- 
tron chains there are important deviations with respect 
to this rule, as for instance in a-decay from odd mass nu- 
clei, and this feature is strongly connected with nuclear 
structure details. Let us mention that different forms of 
the Viola-Seaborg relation were used in Refs. @, 0). 

The Viola-Seaborg rule can be generalized for heavy- 
cluster decays [8j, as it is shown in figure [2l Here the an- 
gular momenta carried by emitted fragments are zero. In 
Ref. 0] it was proposed the following generalized Viola- 
Seaborg rule for the heavy cluster emission 



log 10 T = ai 



ZdZc 



a-iZ^Zc 



C-2 



(2.11) 



with the following set of parameters 

oi = 1.517 , a 2 = 0.053 , b 2 = -92.911 , 
c 2 = {even), = 1.402 (odd) , 

where c 2 is the blocking parameter for odd-mass nuclei. 
Thus, from Eq. (|2.8p the reduced width in this case 
should depend upon the product between daugher and 
cluster charges. 

An interesting feature can be seen by plotting in fig- 
ure [3] the logarithm of the reduced half life (|2.8[) ver- 
sus the Coulomb parameter for various proton emitters 
[ToL [ll[ . In this case most of emitters have non- vanishing 
angular momentum. The data are rougly divided into 
two regions, corresponding to Z < 68 (open circles) and 
Z > 68 (dark circles). This pattern can be asimilitated 
with a generalized Viola-Seaborg rule for the two groups 
of charge numbers. 

The situation with binary cold fission is quite different. 
Here there is not a simple linear dependence between the 
half-life for a given isotopic partition and the Coulomb 
parameter, due to the fact that mass asymmetry changes 
during the scission process. 



3. A SIMPLE MODEL FOR THE REDUCED 
WIDTH 

The decay processes can be schematically described by 
the following cluster-daugher spherical potential 

V(r) = fa / (r ~ r ° )2 +v , r<r B 
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FIG. 3: Logarithm of reduced half lives i2.8\) for proton emit- 
ters versus Coulomb parameter (2.30 . Open circles denote 
emitters with Z < 68, while dark circles emitters with Z > 68. 
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where ro = 1.2 Ap is the surface radius of the daugher 
nucleus. Indeed, microscopic calculations have shown 
that the preformation factor of the a-particle has a Gaus- 
sian shape, centered on the nuclear surface [l2[. More- 
over, the spherical component of the prefomation am- 
plitude gives more than 90% contribution in the a-decay 
width of deformed nuclei. Notice that the radial equation 
of the shifted harmonic oscillator (ho) potential is similar 
with the equation of the one-dimensional oscillator, but 
having approximate eigenvalues given by 



E„ 



i\ h 2 l{l + \) 



(3.2) 



By considering Q-value as the first eigenstate in the 
shifted ho well Q — vq = together with the conti- 

nuity condition at the top of the barrier rs, one obtains 
the following relation 

(3(r B -r ) 2 , . 1 

hw = V frag (r B ) + -JUo , (3.3) 

where it was introduced the so called fragmentation (or 
driving) potential, as the difference between the top of 
the Coulomb barrier and Q-value 



V frag {r B ) = V c {r B ) - 



(3.4) 



V c {r) 



r > r B 



(3.1) 



The second component in Eq. (|2.8[) contains the 
logarithm of the Coulomb-Hankel function inside the 
Coulomb barrier which, according to Eq. (|2.6p . is pro- 
portional with the Colomb parameter \- The third part 
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contains the reduced width squared which, according to 
Eq. (|2.2[) . is proportional with the modulus of the inter- 
nal wave function squared. For a shifted ho well one has 
for the ground state 



\f^(r)\ 2 =Ala-^- r °r . 



(3.5) 



By using the notation 7 = 70 and Eq. ()3.3() one obtains 
the following relation 



4. DECAY RULE FOR REDUCED WIDTHS 

Most of experimental data refer to the a-decay. There- 
fore there were analyzed reduced widths in a-decays con- 
necting ground states of even-even nuclei. In figure[?]it is 
plotted the logarithm of the experimental reduced width 
squared, by using the above relation versus the fragmen- 
tation potential. The data are divided into five regions 
of even-even a emitters as follows 



2 log w e z h 2 A^ 

logi 7 (rs) = j—Vf rag(rB) +l0g W 7j^^ (3-6) 
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FIG. 4: The logarithm of the a-decay reduced width sq 
versus the fragmentation potential f or regions of the nu- 
clear chart described by fJT/p. 



In this way one obtains that indeed the logarithm of 
the half life is of the Viola-Seaborg type 

lo gio T = c i( r s)x + c 2 Vf ra g(r B ) + c 3 (r B , A%) , (3.7) 

because the fragmentation potential contains the product 
ZryZc- Its coefficient depends upon the touching radius, 
but this radius has a very small variation along various 
isotopic chains. Notice that the slope in Eq. (|3.6[) has 
a negative value and it is connected with the shape of 
the interaction potential (ho energy hui), while the free 
term gives information about the amplitude of the cluster 
wave function. 

Our calculation has shown that the linear relation p.6p 
but with different coefficients, remains valid in the most 
general case of the double folding plus repulsive interac- 
tion between fragments, used in Refs. [HI, EH- 



1) Z<82, 50<A< 82 Fig. H(e), stars ; (4.1) 

2) Z< 82, 82<A<126 Fig. H (a), crosses ; 

3) Z> 82, 82<A<126 FigMQ>), circles ; 

4) Z > 82, 126 < N < 152 Fig.\Q(c), squares ; 

5) Z > 82, A>152 FigM(d), triangles . 

In calculations it was used the value of the touching 
radius, i.e. 



tb 



1.2(^ 3 + ^H , 



l/3x 



(4.2) 



where the approximation Kn w 8w in Eq. (|2.4[) is ful- 
filled with 90% accuracy for the most deformed nuclei 
[1]. Notice that regions 1-4 contain rather long isotopic 
chains, while in the last region 5 one has not more than 
two isotopes/chain. This is the reason why, except for the 
last region 5, the reduced width decreases with respect 
to the fragmentation potential, according to the theoret- 
ical prediction given by (|3 . 6|) . Notice that in the regions 
1 and 4, above 100 Sn and 208 Pb double magic emitters, 
respectivelly, the ho parameter of the a-daughter poten- 
tial is larger than for regions 2 and 3, corresponding to 
charge numbers around the double magic nucleus 208 Pb. 
On the other hand, one obtains the largest amplitudes of 
the cluster wave function in the regions 2 and 3. 

An interesting observation can be made from figure [5l 
where it is plotted the difference log 10 T — c-xVfragiTB) — 
c 3 (rB,AQ) versus the Coulomb parameter \, by using 
the same five symbols for the above described regions. 
Amazingly enough there were obtained three lines, cor- 
responding to different amplitudes of the cluster wave 
function. The regions 1 and 4, corresponding to emitters 
above double magic nuclei 50 Sn and 208 Pb, respectivelly, 
have practically the same internal amplitudes Aq. The 
same is true for the regions 3 and 5. 

The linear dependence of logio^f 2 versus the fragmenta- 
tion potential (|3.6p remains valid for any kind of cluster 
emission. This fact is nicely confirmed by heavy clus- 
ter emission processes in figure M (a) , where it is plotted 
the dependence between the corresponding experimental 
values for the same decays in figured] Here it is also plot- 
ted a similar dependence for a-decays corresponding to 
the same heavy cluster emitters. The straight line is the 
linear fit for cluster emission processes, except a-decays 



log W "f 2 = -0.586(V C - Q) + 15.399 



(4.3) 
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FIG. 5: The difference log 10 T — C2Vf rag {rB) — C3 versus the 
Coulomb parameter \ f or fi ve different regions described by 
\4- The straight lines are the corresponding linear fits. 



FIG. 6: (a) The logarithm of the reduced width squared versus 
the fragmentation potential \3.4\ ). Different symbols corre- 
spond to cluster decays in figure The straight line is the 
linear fit for cluster emission processes, except a-decay. 
(b) Cluster mass number versus the fragmentation potential. 



The above value of the slope —logioe 2 /fko in Eq. ()3.6|) 
leads to ftu « 1.5 MeV, with the same order of magnitude 
as in the a-decay case. The relative large scattering of 
experimental data around the straight line in figure[6]can 
be explained by the simplicity of the used cluster-core 
potential (|3.1[) . 

Let us mention that a relation expressing the spectro- 
scopic factor (proportional with the reduced width) for 
cluster emission processes was derived in Ref. [l5| 

S = S^" 1 )/ 3 , (4.4) 

where Ac is the mass of the emitted light cluster and 
S a ~ 10~ 2 . As can be seen from figure [5] (b) between 
Ac and Vf rag there exists a rather good linear depen- 
dence and therefore the above scalling law can be easily 
understood in terms of the fragmentation potential. 

Concer ning the reduced widths of proton emitters in 
Refs. [Til flq] it was pointed out the correlation between 
the reduced width and the quadrupole deformation. This 
fact can be seen in figure[jj(a), where the region with Z < 
68 corresponds to f3 > 0.1 (open circles), while the other 
one with Z > 68 to j3 < 0.1 (dark circles). The two linear 
fits have obviously different slopes. This dependence is 
induced by the propagator matrix Kiy (/3, r) in Eq. (|2.4|) . 
Notice that the two dark circles with the smallest reduced 
widths correspond to the heaviest emitters with Z > 80. 

At the same time one sees from figure [Jj (b) that the 
same data are clustered into two regions, which can be 
directly related with the fragmentation potential (|3.4p . 
Here the two linear fits in terms of the fragmentation po- 



tential, corresponding to the two regions of charge num- 
bers, have roughly the same slopes, but different values 
in origin. Thus, the two diferrent lines in figure [3] can 
be directly connected with similar lines in figure [7J (b) . 
They correspond to different orders of magnitude of the 
fragmentation potential, giving different orders to wave 
functions and therefore to reduced widths. 

A special case is given by the two-proton emission. 
This process was predicted log time ago [171 ]. but only 
few such emitters were recently detected until now. Let 
us mention that the most recent general treatment of the 
two-proton emission process, assuming a three-body dy- 
namics, is given in Refs. [l8j (and the references therein). 

The experimental half-lives versus Coulomb parameter 
are given in figure[5](a) by triangles. The emitter charges 
are also pointed out. Here it is assumed a simplified ver- 
sion, where the light emitted cluster is supposed to be 
a di-proton with the charge Zc = 2. In this case one 
can use the factorisation of the decay width (|2.ip . One 
sees that half lives (triangles) follow the general trend 
(dashed line) of the usual proton emitters, given in the 
same figure by the symbols in figure [3] In figure [8] (b) 
it is plotted the logarithm of the reduced width squared 
versus the fragmentation potential by triangles. One in- 
deed observes that the slope of the fitting dashed line has 
a negative value — logioe 2 /Hui, but with a much larger 
value in comparison with proton emitters, given by the 
same symbols as in figure [7J (b) . The three lines in this 
figure, corresponding to proton and two-proton radioac- 



6 




9 9.5 10 10.5 11 11.5 12 12.5 13 13.5 14 2 4 6 8 10 12 14 



FIG. 7: (a) The logarithm of the reduced width squared versus 
the quadrupole deformation. By open circles are given emit- 
ters with Z < 68, while by dark circles those with Z > 68 
for proton emission. The two regression lines fit the corre- 
sponding data, (b) The logarithm of the reduced width squared 
versus the fragmentation potential {3.$ - The symbols are the 
same as in (a). 



FIG. 8: (a) The logarithm of the half life versus Coulomb 
parameter for two-proton emitters (triangles). By circles are 
given data for proton emitters in figure^ (b) The logarithm of 
the reduced width squared versus the fragmentation potential 
\3. 4\ ) for two-proton emitters (triangles) . The same quantity 
is given by circles for proton emitters in figure^(b). The solid 
lines fit proton data, while the dashed lines fit two-proton data. 



tivity, respectively, are given by 

log wl 2 = -0.283(V C - Q) + 1-329 , Z < 68 
log wl 2 = -0.365(V C - Q) + 3.440 , Z > 68 
logwl 2 = -2.075(Vb -Q) + 4.403 . (4.5) 

The ho energy is Huj w 1.5 MeV for proton emission (i.e. 
the same order as for heavy cluster radioactivity and a- 
decay) and fuv ~ 0.2 MeV for two-proton emission, by 
considering the di-proton approximation. This small ho 
energy is connected with the use of Eq. (|4.2p in estimat- 
ing the spatial extension of the di-proton system Rq and 
therefore the fragmentation potential. In reality the di- 
proton is not a bound system and it changes its size dur- 
ing the barrier penetration. Our microscopic estimate, 
by using a pairing residual interaction between emitted 
protons from 45 Fe, evidenced that the size of the wave 
packet increases in the relative coordinate by 1 fm over 
a distance of 1 fm in the region of the nuclear surface. 
Actually the reduced width defined by (|2.ip can be gen- 
eralized by using the other extreme scenario, given by a 
sequential emission, where this relation is integrated over 
all possible energies of emitted particles [l||. 

An interesting observation concerns the amplitude Aq 
in Eq. (|3.6p . given by the value of fitting lines in ori- 
gin. It has similar values for both two-proton and proton 
emitters with Z > 68. 



Now let us analyze a-decay processes to excited low ly- 
ing 2 + states. There were considered more than 70 decays 
of even-even rotational, vibrational and transitional nu- 
clei [ill EH • The hindrance factor (HF) is defined as the 
ratio of reduced widths squared connecting the ground 
states and ground to excited states with the angular mo- 
mentum / = 2, i.e. 

HF=^% . (4.6) 

72 

Thus, by using (|3.6[) . the logarithm of the HF becomes 
proportional with the excitation energy of the daughter 
nucleus 

log w HF = — E y 2 + l°9io^2 • ( 4J ) 

It is worth mentioning that this relation is equivalent 
with the Boltzman distribution for the reduced width to 
the excited state 72. In Refs. [l^, [2l[ such a dependence 
was postulated in order to describe HF's. 

In figure [5] it is plotted the logarithm of the HF versus 
the excitation energy for rotational nuclei with E^ x ^ < 
0.1 MeV, by using the same notations given by Eq. (|4.ip . 
As a rule the HF's have small values and therefore the 
wave functions have similar amplitudes io ~ ^2- Notice 
that the region 4 in figure [HI (b), with Z > 82, 126 < 
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FIG. 9: The logarithm of the hindrance factor versus the ex- 
citation energy of the daughter nucleus for rotational nuclei. 
The symbols and numbers correspond to the regions given by 
Eq. Ojp. 



of magnitude, corresponding to an ho energy hu> ~ 1.5 
MeV, for all decay processes, except the di-proton emis- 
sion, where the factorisation (|2.ip is not anymore valid. 
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FIG. 10: The same as in figure but for transitional ana 
vibrational nuclei. 



N < 126, contains most of rotational emitters (33 out 
of 41). Moreover, our analysis has shown that here the 
HF has an almost constant value along various isotopic 
chains, due to the fact that the energy range is very short 
(about 100 keV). 

In figure [TO] it is given the same quantity, but for tran- 
sitional and vibrational nuclei, with > 0.1 MeV. 
The situation looks here to be different and more com- 
plex, with respect to rotational nuclei. The best example 
is given by the same region 4 in figure [TD] (c) , where the 
slope has a positive value, as predicted by Eq. (|4.7|) . No- 
tice that this region contains almost half of the analyzed 
vibrational emitters (14 out of 31). 

5. CONCLUSIONS 

Cold emission processes in terms of the well known fac- 
torisation of the decay width between the penetrability 
and reduced width squared, were analyzed. Based on a 
simple model of the two-body dynamics, namely a shifted 
harmonic oscillator potential surounded by the Coulomb 
interaction, it was derived an universal analytical relation 
expressing the logarithm of the reduced width squared as 
a linear function in terms of the fragmentation potential, 
defined as the difference between the Coulomb barrier 
and the Q- value. Notice that the slope has the same order 



This rule is a consequence of the fact that the logarithm 
of the wave function squared is proportional with the dif- 
ference between the height of the Coulomb potential at 
a given radius and the energy of the system. It is ful- 
filled with a reasonable accuracy by experimental data, 
describing transitions between ground states as well as 
for a-transitions to excited states in the most relevant 
region with Z > 82, 126 < N < 152. As a particu- 
lar case, the two regions in the dependence of proton 
emitters half-lives, corrected by the centrifugal barrier, 
versus the Coulomb parameter are directly related with 
the corresponding regions of the fragmentation potential. 
Thus, the clustering character of the wave function inside 
the Coulomb barrier (i.e. the fragments are already pre- 
formed here) is evidenced by this linear rule in terms of 
the fragmentation potential. On the other hand the well 
known scalling relation, connecting the spectroscopic fac- 
tor with the mass of the emitted cluster, can be nicely 
explained in terms of the fragmentation potential. 
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